Spin amplification, reading, and writing in transport through anisotropic magnetic molecules 
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Inelastic transport through a single magnetic molecule weakly coupled to metallic leads is studied theoreti- 
cally. We consider dynamical processes that are relevant for writing, storing, and reading spin information in 
molecular memory devices. Magnetic anisotropy is found to be crucial for slow spin relaxation. In the presence 
of anisotropy we find giant spin amplification: The spin accumulated in the leads if a bias voltage is applied to 
a molecule prepared in a spin-polarized state can be made exponentially large in a characteristic energy divided 
by temperature. For one ferromagnetic and one paramagnetic lead the molecular spin can be reversed by apply- 
ing a bias voltage even in the absence of a magnetic field. We propose schemes for reading and writing spin 
information based on our findings. 

PACS numbers: 73.63.-b, 75.50.Xx, 85.65.+h, 73.23.Hk 



I. INTRODUCTION 

Recently, a number of fascinating experiments have studied 
transport through magnetic molecules. The current through a 
molecule weakly coupled to metallic leads shows steps as a 
function of bias voltage when additional molecular transitions 
become energetically available. The differential conductance 
dl/dV then shows peaks. Zeeman splitting of the energy lev- 
els and of the peaks in dl/dV has been observed for mag- 
netic moleculesJi^ Magnetic anisotropy partially lifts the de- 
generacy of molecular levels, leading to fine structure of the 
peaks even for vanishing magnetic field, which has been ob- 
served for Mni2 complexes?^ Some of the fine structure peaks 
show negative differential conductance (NDC)i^ The Kondo 
effect has also been found for single moleculesil^ The Zee- 
man splitting"* as well as spin blockade and NDCA- in mag- 
netic molecules have been investigated theoretically. Trans- 
port through magnetic molecules has also been studied in the 
STM geometry, where the coupling to the substrate is much 
stronger.^ 

Magnetic molecules are interesting for applications that 
combine spintronicspiS i.e., the idea to use the electron spin 
in electronic devices, with molecular electronics, i.e., the uti- 
lization of molecules as electronic componentsi^iiS We study 
processes that are relevant for the utilization of magnetic 
molecules as memory cells, which has been discussed at least 
since the early 1990si2iii4i2ii^ In this context it is crucial that 
the stored information persists over sufficiently long times and 
that one can read out and write the information. The read 
out process would naturally utilize the interaction of electrons 
tunneling through the molecule with its local spin. The works 
mentioned above lay the ground for a study of this aspect. 
On the other hand, the other two aspects — ^persistence and 
writing — have received little attention. 

In the present paper, we consider current-induced spin re- 
laxation and propose schemes for reading and writing the 
molecular spin electronically. We show that the spin trans- 
mitted from one lead to the other before the molecular spin 
relaxes can be made very large. This interesting effect of spin 
amplification forms the basis of our proposed read-out mech- 



anism. Magnetic anisotropy is a key ingredient for this effect 
and for molecular memory devices in general. Since in the ab- 
sence of a magnetic field the spin polarization vanishes in the 
steady state, we expect any spin polarization to decay in time. 
For memory applications this relaxation should be slow. That 
this can be accomplished by an easy-axis anisotropy, which 
introduces an energy barrier for a change of spin direction?^ 
has been shown for Fes complexes by Barra et alM Typical 
molecules we have in mind are planar complexes of magnetic 
ions, e.g., porphyrin complexes Uke heme. 



II. THEORY 

We consider a magnetic molecule with one orbital coupled 
to a local spin S, described by the Hamiltonian 



U 



Hmoi = (e - eVg) n+ — h{n - 1) - J s • S 



(1) 



Here, n 



- J 



jCf + c|c| is the number operator {c\ creates an 

electron in the orbital) and s = X^o-o-' clj{crcrcr' /'^)ca' is the 
corresponding spin operator, e is the single-electron energy in 
the orbital, which can be shifted by the gate voltage Vg, U is 
the Coulomb repulsion between two electrons in the orbital, 
and J is the exchange interaction between an electron in the 
orbital and the local spin of length S >1. Easy-axis magnetic 
anisotropy is assumed {K^ > 0) for the local spin. We define 
the total spin operator Stot = s + S. 

The Hamiltonian in Eq. Q is the simplest one with an 
anisotropy-induced energy barrier and already exhibits the es- 
sential physics of reading, storing, and writing molecular-spin 
information, as we show in this paper. Our model does not 
allow for magnetic tunneling through the energy barrier in the 
absence of electron tunneling involving the leads. This effect 
is included in Ref. '5 but is assumed to be weak there. Our 
model is valid as long as the typical rate of spontaneous mag- 
netic tunneling is small compared to typical electronic tunnel- 
ing rates. 
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The molecule is weakly coupled to two leads L (left) and R 
(right) with Hamiltonians Ha = I]ka '^"ktr a^ka^akcr where 
"^QkiT creates an electron in lead a = L, R. The hybridization 
of the molecular orbital with the leads is described by Hamil- 
tonians iJhyb^a = Ek^rl^ 

a ^akcr^^ li.c). In the following 
we consider the case of symmetric contacts, = t-R,, to avoid 
unnecessary complications that have nothing to do with the 
physics we are interested in. However, it should be noted that 
asymmetry of the contacts, which is expected to be present, 
can have a significant effect on transport properties, 

The unperturbed Hamiltonian i^moi is easily diagonal- 
ized. The eigenstates fall into sectors with n = 0,1,2 
electrons. Since [S'fo,. , i?moi] = the eigenvalue m of 
Sfot is another good quantum number For n = 0, 2 we 
obtain 2(25 + 1) eigenstates {|0, m), |2,m)} with energies 
e(0, m) = -K2m^ - Bm and e(2, to) = 2(e - eVg) + [/ - 
K2'm? — Bm, respectively. For ri = 1 we find the eigenstates, 
in terms of spin states of the electron and the local spin, 



J) TO 



|l,77l)=^ =T — 



2VA£; 
JVS'(5 + 1)-to2 + 1/4 



2^fKEyj2/\E =F (2A'2 - J)m 
with energies 

J 



i>|m + l/2) 
T)|to-1/2), (2) 



e=^(l,TO) = e~eVg-Bm+--K2\m'' + -'j±/^E{'m) (3) 

where l\E{m) = [K2{K2 - J)m^ + {J/Af(2S + ifY'"^. 
There are two orthogonal states 1 1 , to) ^ for — 5 + 1/2 < to < 
5 - 1/2. For the fully polarized states |l,-5 - 1/2) and 
1 1,5 + 1/2) the upper (lower) sign applies if A'2 — J/2 is 
positive (negative). 

The hybridization is treated as a perturbation. The deriva- 
tion follows Refs. i4»16b17.181 and is not repeated here. We 
obtain coupled rate equations 



pi 



(4) 



for the probabilities P* of the molecular many -body states |i) . 
The steady state Pq is the solution for = 0. 

The transition rates Ri^j are written as a sum over the two 
leads and the two spin directions, Ri^j = J2cra ^-i-^j' ^^'■^ 



To 



■eSaV/2)\C, 



a |2 



eSaVl2)\Clf\ 



(5) 



where sl.r = ±1, I/tq = 27r |tp Dvhc/Tt- is the typical tran- 
sition rate in terms of the density of states D (for one spin 
direction) of the leads and their unit-cell volume Uuc, V is the 
bias voltage, e,; is the energy of state |i), f{x) is the Fermi 
function, and C^"^- = {i\ca\j)- 

The occupation and spin in a state described by probabili- 
ties are given by n = riiP^ and to = miP\ re- 
spectively. Here, Ui (nii) is the occupation (polarization) in 



state \i). The current through lead a is 



){R]Z. + RjZ.) P\ 



(6) 



where we count currents from left to right as positive. The 
steady -state currents (/") through the two leads are equal 
since the occupation is constant in time. The current scale 
is g/tq, where the rate I/tq comes from the transition rates. 
Comparing this to typical experimental currents of 0.1 nAr*iS 
we obtain tq ~ 1.6 ns. The spin current 



R 



la 



(7) 



is generally nonzero if the molecule is prepared in a spin- 
polarized state. The molecule relaxes exponentially towards 
the steady state, which for P = is nonmagnetic. We can 
thus define the total transmitted spin 



A5"[P] 



dti:it,p), 



(8) 



which depends on the probabilities P ~ (P^, . . .) at t = 0. 

To evaluate Eq. (|8} we require the time dependence of P. 
We rewrite the rate equations in matrix form, P — AP, 
where A has the components Aij = Rj-,i for i j and 



R, 



The solution is P{t) 



r,At 



P(0), 



which allows to calculate occupation, polarization, and (spin) 
current as functions of time for given initial conditions. 



in. RESULTS 

As noted above, the anisotropy partially lifts the degeneracy 
of molecular levels. The resulting fine structure of the peaks 
in dl/dV is shown in Fig.fQa) for vanishing anisotropy and 
in Fig.^b) for K2 = 0.04. In the latter case one of the peaks 
exhibits NDC. To elucidate this effect we plot in Fig. [Qc) 
the current, differential conductance, and relevant occupation 
probabilities at constant gate voltage at a lower temperature. 

At low bias voltage only the degenerate ground states 
|l,5/2) and |1,— 5/2) are occupied and no current is flow- 
ing. For increasing bias first the transitions to |0, ±2) be- 
come possible, cf. the energy-level scheme in the inset in 
Fig- Wic). The current then increases to a plateau. On this 
plateau the four states have equal probability and Eq. (jSJl 
yields |(/'")| = § ±5/2)^10 ±2)- Next, the transitions 
from |0, ±2) to |l,±3/2)~ become active and the system 
reaches a new plateau, denoted by an arrow in Fig.^c), with 



— o (^|i,±5/2>-»|0,±2) + -fi^|l,±3/2)-^|0,±2))- (9) 



Since 



nL 

|l,±3/2) 



particular parameters 
< 1/2 due to the 



for these 

♦ |0.±2)/-^l,±5/2)^|0,±2> 

coefficients in Eq. (|2}, the current is smaller than on the 
first plateau, leading to NDC in the crossover region. The 
transitions between plateaus are rounded due to the nonzero 
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FIG. 1: (Color online) dl /dV as a function of gate voltage Vg and 
bias voltage V . Yellow (bright) colors denote positive values of the 
differential conductance, dark red (dark gray) colors denote nega- 
tive values. In (a) results for an isotropic molecule are shown, the 
parameters are T = 0.002, 5 = 2, e = 0.2, [/ = 1, J = 0.1, 
K2 = B — 0. Energies are here and in the following given in elec- 
tron volts. Coulomb blockade is found to the left and right of the V- 
shaped region, whereas within this region a nonzero steady-state cur- 
rent is flowing. The satelite line results from the exchange splitting 
of energy levels for n = 1. In (b) an anisotropy of K2 =0.04 has 
been assumed, leading to a complex splitting of the dl /dV peaks. 
We find one peak with NDC (arrow), (c) Occupation, steady-state 
current in units of e/ro, dl /dV in units of e/ro per volt, and prob- 
abilities of various states (dashed) for a cut through the NDC region 
in (b) at the lower temperature T — 0.001. The arrow indicates the 
plateau with reduced current. Inset: Molecular energy levels of states 
with n — (black bars) and n — 1 (red/gray bars with circles) and 
magnetic quantum numbers m. 



temperature. This mechanism is different from the one for 
Mni2 complexes^ — in our case there is no blocking state 
in which the molecule becomes trapped due to suppressed 
outgoing transition rates. On the contrary, the probabilities 
P' are equal for all accessable states. 

We now turn to the relaxation of the molecular spin. For a 
molecule prepared in state |0, 2) at time t = 0, Fig.|2ja) shows 
the time dependence of occupation and spin polarization. The 
occupation approaches the constant value (n) on the timescale 
To, the typical time for a single tunneling event. 

The spin polarization shows a quite different behavior with 
two distinct timescales. Initially, m{t) approaches a "quasi- 
steady" state on the timescale tq, which in this case has higher 
polarization since the state 1 1,5/2) has significant weight. 
Then, m(t) decays to zero more slowly. This decay is slow 
because the molecule must pass through several intermedi- 
ate states to reach a state with opposite spin polarization, es- 
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FIG. 2: (Color online) (a) Occupation n{t) and spin polarization 
m{t) for bias V = 0.2 V, and T = 0.002, S = 2, e = 0.12, 
(7 = 1, J = 0.1, and B = for various anisotropics A"2. In this and 
the following figures we absorp the gate voltage into e. The molecule 
is in state |0, 2) at time t = 0. (b) Molecular energy levels for states 
with n = (black bars) and n — 1 (red/gray bars with circles) for 
the parameters from (a) and K2 = 0.04. Active transitions are de- 
noted by dashed lines. The crosses denote the transitions that next 
become thermally suppressed for larger K2 or smaller bias V. 



sentially performing a one-dimensional random waUc, as indi- 
cated in the level scheme, Fig.|2jb). Figure|2ja) shows that the 
spin relaxation initially becomes slightly /aifer for increasing 
K2- This is due to the change of matrix elements C[j with 
K2- For anisotropics K2 ^ 0.04 the decay becomes very 
slow since two of the transitions needed to reverse the spin, 
denoted by the crosses in Fig.|2jb), become higher in energy 
than eV/2 and are thus forbidden for T ^ and thermally 
activated for T > 0. 

Importantly, in the regime of thermally activated (slow) 
spin relaxation a sizeable steady-state current is flowing since 
transitions between the states |0,±2) and |l,±5/2) are still 
possible, see Fig.|2jb). As noted above, this should lead to a 
nonzero total transmitted spin, Eq. (|8j. This is indeed found 
in Fig. 13 Varying the bias voltage we observe four regimes: 

(i) For small bias we are in the Coulomb blockade regime 
with very small steady-state current. However, we find that 
the total transmitted spin is exponentially large in a charac- 
teristic energy barrier AE over temperature. In our example 
AE is the difference between the energy e^(l, 3/2) — e(0, 2) 
necessary for a spin-down electron to tunnel in and the avail- 
able energy eV/2 of an incoming electron, which is smaller. 
While all transitions are thermally activated, the ones neces- 
sary to relax the spin have a much higher energy than the 
transition between |0, 2) and |l,5/2), which dominates the 
current. However, the transmission takes exponentially long 
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FIG. 3: (Color online) Total transmitted spin through the left and 
right leads for a magnetic molecule prepared in the n = maximum 
spin state |0, 2), as functions of bias voltage. The parameters are as in 
Fig.|2|with A'2 ~ 0.04. The steady-state current is also shown. The 
inset shows the spin transmitted through the left lead on a logarithmic 
scale as a function of bias voltage for three different temperaturesi— 



since this transition is also thermally activated. 

(ii) For larger bias we find the most interesting regime. The 
spin relaxation rate is still small while the current is large. The 
bias is too small to overcome the energy barrier between spin 
up and down, see Fig.|2jb). If the systems starts in state |0, 2), 
the only transition with large rate is to 1 1 , 5/2) ; the transitions 
to |1, 3/2)^ are thermally suppressed in this regime [note that 
they are higher in energy in Fig.|2jb)]. Thus a spin-up electron 
has to enter since the process increases the total spin from 2 
to 5/2. This electron can leave the molecule through the other 
lead, returning it to the intial state. On the other hand, a sec- 
ond electron cannot enter the molecule since this would cost 
a high Coulomb energy U. Therefore, the current if fully spin 
polarized until a transition to |1, 3/2)^ occurs. This thermally 
suppressed transition happens with a small rate proportional to 
exp(-A£;/T), where AE = £-(1, 3/2) - e(0, 2) - eV/2, as 
above. Consequently, the current is spin-polarized for an ex- 
ponentially long (in AE/T) time leading to an exponentially 
large transmitted spin |A5"| [the negative sign in Fig.|3lcan 
be understood by considering the transition rates, Eq. (|5j, in 
detail]. The exponential dependence of the transmitted spin 
on AE/T is clearly seen in the inset in Fig.|3] The average 
time = l/|/"(< = 0)| for one unit of spin to be trans- 
mitted is of the order of tqi^ If is short compared to the 
spin relaxation time in the leads large opposite magnetizations 
will be accumulated in the leads, similarly to the effect of spin 
injection studied intensively in recent yeaisM This effect is 
induced by the breaking of spin symmetry at t = only 
through the polarization of a single quantum spin, and can 
thus be described as giant spin amplification. It is a promis- 
ing method to read out the spin information. The magnetiza- 
tion in the leads could be detected with a pickup coil or by the 
magneto-optical Kerr effecti2i2i A strong amplification mech- 
anism could also facilitate the reliable transfer of spin between 
individual molecules in a device. 



initial state: 10, -2> 




time t/Xg 

FIG. 4: (Color online) Probabilities of molecular states after initial 
preparation in state |0, — 2) as functions of time at the bias voltage 
V — 0.207 V. The probabilities of states with extremal spin polar- 
ization and of the representative intermediate state jO, 0) are shown. 
A magnetic field (B = 0.002) leads to spin relaxation towards larger 
spin values. The parameters are T — 0.0002 (lower than before), 
K2 = 0.04 and otherwise as in Fig.|3a). Inset: Probabilities after 
the bias voltage is suddenly switched off in the steady state reached 
for V = 0.207 V. 



(iii) Further increase of the bias leads to a regime where the 
transmitted spin is nonzero but not exponentially enhanced. 
Here, both spin relaxation rate and current are large. The 
bias is large enough to overcome the energy barrier. Figure 
|3] shows that the spin transmitted through the left lead is (in 
absolute value) larger by two than through the right lead. The 
reason is that the electrons flow from right to left and that the 
spin of the initial state has to leave the molecule for t — > 00. 

(iv) At large bias all transitions between states with n = 
0, 1 are possible. The spin transmitted through the right, 
incoming lead, AS^, essentially vanishes. The initial spin 
leaves the molecule through the left lead, leading to AS^ = 
—2. This is the nonmagnetic regime. 

The pattern seen in Fig. |3] is robust under change of 
anisotropy A'2. For K2 = the regime of giant spin am- 
plification is absent but for any K2 > it exists at sufficiently 
low temperatures since the only requirement is that the energy 
of state 1 1 , 5/2) is larger than the one of 1 1 , 3/2) ~ . 

Finally, for molecular memory application one also needs 
to write the information, i.e., to switch the molecule to a pre- 
determined state. An obvious idea is to apply a magnetic field 
to a molecule attached to the nonmagnetic leads considered so 
far However, applying a field at zero bias does not work since 
all transitions remain thermally suppressed so that relaxation 
to the spin-polarized steady state is exponentially slow. 

We find that reliable switching requires a two-step scheme: 
First one applies a magnetic field, which tilts the energy lev- 
els in Fig. |3b), and a bias voltage that is just large enough 
to allow transitions in the desired direction but not in the op- 
posite one at low temperatures. Since the Zeeman energy B 
is small this requires fine tuning of V on the scale B/e and 
cooling to T <^ B. Figure |3 shows the change of proba- 
bilities for all states with n = 0, 1 if one starts with |0, —2) 
and applies a positive B field to switch the molecule to spin 
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FIG. 5: (Color online) Probabilities of molecular states after initial 
preparation in state (a) |0, —2) and (b) |0, 2) as functions of time 
at the bias voltage (a) 1^ = 0.7 V and (b) -0.7 V. The right lead is 
ferromagnetic with a ratio of densities of states of j — 0.001. 
There is no applied magnetic field. The other parameters are as in 
Fig.|4|except for the much higher temperature T = 0.02. 



up. The Zeeman energy is chosen as 2 mcV. The molecule 
crosses over from |0,— 2) to a steady state essentially con- 
sisting of |0, 2) and |1, 5/2) on a timescale of the same order 
and of the same origin as the spin relaxation times, see Fig.|2] 
Note that in a molecular circuit one could apply a magnetic 
field to many molecules and address a specific one with the 
bias voltage. In the second step the bias is switched off. The 
inset of Fig.0]shows that the system then relaxes towards the 
target state |0, 2) on the timescale tq. 

This writing scheme is problematic because it requires 
to reverse a large magnetic field between switching events, 
which is a very slow process. It also requires very low temper- 
atures. We can overcome these problems by using one ferro- 
magnetic and one nonmagnetic lead. The ferromagnetic lead 
R is modelled by different densities of states for a =t, j 
electrons. We set jD^ = 0.001 so that the lead is essen- 
tially a half-metallic ferromagnet such as NiMnSb and Cr02. 
For less complete spin polarization in the lead the reliability 
of the switching degrades. We note that spin injection from 
a ferromagnet is not trivial. However, in recent years signif- 
icant progress has been maden^^i Nevertheless, the most diffi- 
cult task in implementing this writing scheme probably is to 
achieve a high degree of spin polarization of tunneling elec- 
trons. 

Figure |5] shows that one can switch the spin polarization 
in hoth directions on the timescale of typical spin relaxation 
times by applying a bias voltage in vanishing magnetic field. 
For > 0, Fig. ISja), the electrons flow from right to left. 
Since is small, nearly all electrons have spin up. Due to 



exchange scattering between electron spin and local spin, the 
latter is increased. Figure ^a) shows that all states except 
|0, 2) and |1, 5/2), which have the largest positive magnetiza- 
tion, die out. By switching the bias voltage to = as in 
Fig.0]we can then make the molecule relax rapidly towards 
the unique state |0, 2) (not shown). For V < 0, Fig. |5Jb), 
the electrons flow from left to right. Spin-down electrons are 
essentially trapped on the molecule until they perform a spin 
exchange with the local spin, which decreases the local spin. 
As the result, the state 1 1 , — 5/2) becomes populated at the ex- 
pense of all other states. Importantly, no fine tuning of the bias 
voltage is required. Furthermore, the temperature need not be 
small — Fig.|5]is calculated with T close to room temperature. 
Both properties are desirable for molecular-electronics appli- 
cations. 

If one lead is ferromagnetic, the current is generally spin- 
polarized even in the steady state. Then the total transmitted 
spin, Eq. (|8}, diverges. However, we can define the excess 
transmitted spin in a certain state relative to the steady state, 
AS"[P] = dt[I^{t,P) - (/f)]. This quantity exhibits 
spin amplification as in the case of nonmagnetic leads. 



IV. CONCLUSIONS 

We have studied the inelastic charge and spin transport 
through a magnetically anisotropic molecule weakly cou- 
pled to metallic leads. The three processes crucial for 
molecular memory applications — writing, storing, and read- 
ing information — can be implemented in such a device. The 
information storage is affected by spin relaxation, which can 
be very slow for large easy-axis anisotropy. Also due to the 
anisotropy, application of a bias voltage to a molecule in a 
spin-polarized state can lead to the transfer of a large amount 
of spin or magnetic moment from one lead to the other. This 
transmitted spin increases exponentially for low temperatures. 
We propose that this giant spin amplification could be used 
to read out the spin information. Finally, we also propose 
a scheme to write the information, which does not require a 
magnetic field but uses one ferromagnetic lead. 
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